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ABSTRACT 



Aims. We present a comprehensive study that applies the Fourier transform to a sample of O and early B-type stars (either dwarfs, 
giants, or supergiants) to determine their projected rotational velocities, compare with previous values obtained with other methods, 
and seek for evidence of extra broadening in the spectral lines 

Methods. The Fourier technique, extensively used in the study of cooler stars, has remained only marginally applied for the case of 
early-type stars. The comparison of v sin i values obtained through the ft and fwhm methods shows that the fwhm technique must be 
used with care in the analysis of OB giants and supergiants, and when it is applied to He i lines. Contrarily, the ft method appears to be 
a powerful tool to derive reliable projected rotational velocities, and separate the effect of rotation from other broadening mechanisms 
present in these stars. 

Results. The analysis of the sample of OB stars shows that while dwarfs and giants display a broad range of projected rotational 
velocities, from less than 30 up to 450 kms , supergiants have in general values close to or below 100 fans . The analysis has 
also definitely shown that while the effect of extra broadening is negligible in OB dwarfs, it is clearly present in supergiants. When 
examining the behavior of the projected rotational velocities with the stellar parameters and across the HR diagram, we conclude, in 
agreement with previous researchers, that the rotational velocity should decrease when the stars evolve. On the contrary, macroturbu- 
lence may be constant, resulting therefore in an increasing importance as compared to rotation when the stars evolve. 

Key words. Stars: early-type - Stars: rotation - turbulence 



1. Introduction 

Together with mass and metallicity, stellar rotation (i.e. stel- 
lar angular momentum) determines the formation, structure and 
evolution of a star. Moreover, it is a necessary property for a 
large number of important processes, like the solar dynamo or 
the stellar activity. In the realm of massive stars, rotation is be- 
coming increasingly important. It is thought to be the reason for 
internal mixing mechanisms that would bring CNO processed 
material and transport inner angular momentum to the stellar 
surface (Maeder & Meynet 2000). Coupled with stellar winds, 
it can strongly alter the evolution of the star and the theoreti- 
cal Eddington limit. Equally important, the coupling with stel- 
lar winds imposes a dependence of the rotational velocity itself 
with metallicity during the evolution of high mass stars (Maeder 
& Meynet 2001). Consequently, massive stars in the Magellanic 
Clouds are expected to have larger rotational velocities in aver- 
age, and more CNO contamination in their atmospheres. 

Since present evolutionary stellar models take this rota- 
tional velocity into account, considering it as an important 
parameter in the stellar evolution, accurate determination of 
stellar rotational velocities in massive stars is required, not 
only to constrain these stellar evolutionary models, but also to 
understand the physics of massive stars and their dependence 
with metallicity. 
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Several methods for the determination of the projected rota- 
tional velocity (v sin z) in isolated early-type stars can be found 
in the literature; the most commonly used are: 

1 . The Full-Width-Half-Maximum (fwhm) of the observed line 
profile is measured and related to the v sin i (fwhm method). 

2. The observed spectrum is cross-correlated against a template 
spectrum of low rotational velocity. The Gaussian width of 
the cross-correlation function is directly related to the pro- 
jected rotational velocity of the star if rotation is the domi- 
nant broadening mechanism of the photospheric lines. 

3. The observed line profiles are compared with synthetic ones 
calculated from model atmospheres that are convolved with 
the corresponding broadening functions. Although this is a 
very accurate method if the broadening functions are well 
known, too much computational time is required. Reducing 
that time (e.g., using model grids) reduces also the accuracy. 

4. The Fourier transform of the line profile is calculated. The 
position of zeroes present in this Fourier transform depends 
on the v sin i. 

Being very fast and easy to use, the fwhm method has been ex- 
tensively used for the measurement of projected rotational veloc- 
iti es in O B stars (viz. Sletteb ak [19751 Abt et al. 120021 Herrero et 
al. 119921 Strom et al. 120051 ). The cross-correlation method was 
applied by Penny (1996), Howarth et al. (1997) and Penny et 
al. 020041 ) to iue and stis spectra of O-type stars and early B- 
type supergiants. The studies coming from both methodologies 
(i.e., fwhm and cross-correlation) allow us to have v sin i mea- 
surements of large samples of OB stars with different luminos- 
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ity classes; however, the main limitation of these techniques is 
that they do not allow for the separation of rotational broadening 
from other possible extra broadening mechanisms, usually re- 
ferred to as macroturbulence (although the actual physical mech- 
anism - or mechanisms - remains unknown). In fact, this is a del- 
icate and quite important task, since many works find an absence 
of narrow-lined OB stars. This result has been claimed to be an 
indirect proof of the presence of a non-negligible line broaden- 
ing mechanism in addition to rotation in these stellar objects, 
specially in the case of giants and supergiants (viz. Slettebak 
[1^561 Conti & Ebbets [19771 Howarth et al. [T997t . 

The fitting to broadened synthetic profiles makes it possi- 
ble to distinguish between different broadening mechanisms and 
also to determine other stellar parameters, but as has been com- 
mented before, a large effort is required to use it accurately. 
Furthermore, quite good quality spectra (i.e. high signal-to-noise 
ratio, SNR, and spectral resolution) are required to disentan- 
gle the various broadening mechanisms that may affect the line 
profiles, and the derived values may depend on the considered 
models and the definition of the broadening functions. An early 
attempt to separate the rotational and macroturbulent contribu- 
tions to the line broadening in O- and early B-type stars was per- 
formed by Slettebak (119561 1, however he found that the quality of 
his spectra was not good enough to get any conclusive solution. 

More recently, Ryans et al. (2002) analyzed high quality 
spectra of twelve B-type supergiants. They compare observed 
and synthetically broadened lines using a x 1 technique and find 
for their survey of B-type supergiants that while models domi- 
nated by rotation provide unsatisfactory fits, those where macro- 
turbulence dominates and rotation is negligible (as a broadening 
agent) were acceptable. 

Being also very fast and easy to use, the Fourier method, 
a technique that was pioneered using cool stars, has been only 
marginally applied in the study of OB stars. This technique al- 
lows the v sin i to be easily derived independently of any other 
broadening mechanism affecting the line profile. There are nu- 
merous papers illustrating its use in the determination of v sin i 
in later type stars (viz. Smith |1976al Gray 119801 Royer et al. 
2002a, 2002b, Reiners et al. 2003); however, there is no compre- 
hensive study applying this technique to early-type stars except 
that by Ebbets (119791 ), who analyzed a sample of 16 objects, but 
the limited quality of the spectra allowed him to obtain estima- 
tions of the macroturbulence only for half the objects in the sam- 
ple. Recently, Simon-Diaz et al. ([2006) have applied the Fourier 
method to a small sample of OB stars in the Orion nebula. The 
strength of this method was shown in the study of the main ioniz- 
ing star of this nebula, an 07 V star with a well known rotational 
period of ~ 15 days, that is in conflict with the vsin/ derived 
from the fwhm method, but agrees with the Fourier technique. 

Following the ideas of the present paper and Ryans et 
al. (120021) . Dufton et al. (2006a) have recently analyzed the 
high-resolution spectra for 24 SMC and Galactic B-type super- 
giants to estimate the contributions of both macroturbulence 
and rotation to the broadening of their metal lines. See also the 
analysis of periodically variable B supergiants by Lefever et al. 
(120061 

This paper has three main objectives. First, we extend the 
work presented in Simon-Diaz et al. (120061 1 illustrating the 
strength of the Fourier method in the analysis of early type 
stars; second, we provide new, more accurate values for the ro- 
tational velocities of OB stars and estimate the role played by 
other broadening mechanisms; and third, we extend the work 



by Dufton et al. (2006a) and Lefever et al. (2006) towards early 
spectral types. 

We restored from our archive the spectra of several OB stars, 
obtained in different observing campaigns, for a comprehensive 
determination of the rotational velocities through the Fourier 
analysis method. But firstly, a formal study was done to deter- 
mine the range of applicability of this methodology to OB stars. 

The paper is structured in four main parts as follows: a short 
review of the Fourier methodology, as well as some formal tests 
for the case of massive stars are presented in Sect. [2] The deter- 
mination of the projected rotational velocities of the sample of 
OB type stars is presented in Sect. [3] Finally, the discussion and 
the main conclusions are presented in Sects. [4]and[5] 

2. Testing the Fourier method in early-type stars 

The theoretical base of this methodology was well established 
by Gray ( 119761 see also the latest edition of the book, published 
in 2005). Here, we briefly summarize the main characteristics of 
the method, needed to understand its advantages and limitations, 
and to evaluate the results presented in Sect. [3] We also present 
examples specifically aimed at its application to early-type 
stars. For a more precise description of the methodology and its 
applications to late-type stars, we also refer to Gray's book (see 
above for references). Additional details may be found in the 
pioneering articles by Gray ( 119731 119751 1, Smith ([1976a), and 
Smith & Gray ( 1976b). Some objections to its validity in special 
cases can be found in Mihalas d 19791 ) and Bruning ( 1984). 

In a simple way, an observed stellar line profile may be con- 
sidered to consist of the convolution of an intrinsic profile with 
the rotational, macroturbulent and instrumental profiles. The in- 
trinsic profile includes natural and thermal broadenings, as well 
as other contributions, as those due to microturbulence and Stark 
effect. The Fourier method (ft) for the determination of v sin i 
(c.f. Gray 1 19761 ) is based on the fact that between the previ- 
ous broadening profiles, only the rotational function is expected 
to have zeroes in its Fourier transform (except for large values 
of microturbulence; see below). As firstly described by Carroll 
( 119331 ), the position of these zeroes in frequency space depends 
on the v sin i of the star, so that the frequency of the first zero 
(<T\) is related to the rotational velocity through: 

- v sin i o-i = 0.660 (1) 

c 

Since convolutions transform into products in Fourier space, 
these zeroes will also appear in the total transform function of 
the line profile. Therefore, it is in principle possible to determine, 
independently of any other broadening mechanism, the v sin i of 
a star once the position of the first of those zeroes is identified in 
the Fourier transform of the line profilsfl 

Other mechanisms may contribute with zeroes to the total 
Fourier transform, depending on their profiles in the wavelength 
domain. Note, for example, that microturbulence can also add 
zeroes, when strong lines are considered, due to the boxiness 
caused by it in saturated lines (Gray 1976). However, those ze- 
roes arising from saturation in the Fourier space appear at higher 
frequencies than those for rotation, so there is rarely any confu- 
sion between them, except for those cases with low v sin i. 

1 This argument is only valid under the hypothesis considered by 
the convolution method, which assumes that the effects of the differ- 
ent broadening mechanisms act independently to broaden the emergent 
flux profile (Mihalas [T979l >. 
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2.1. Some formal tests for massive stars 

The properties and advantages of the Fourier method, particu- 
larly to distinguish among physical processes that produce very 
similar profiles in the wavelength domain, have been already dis- 
cussed and illustrated in the works mentioned above. Here we 
present a few applications to the particular case of early-type 
stars. 

For this formal tests we used Fastwind (Santolaya-Rey et 
al. 1 1 9971 Puis et al. 120051 ) theoretical profiles from a model 
with T e ff = 30000 K, log g- 4.0 (an early B dwarf) convolved 
with different broadening profiles. Thermal and microturbulence 
broadenings were included in the line calculation in the usual 
way. 

For the rotational broadenings we considered the usual for- 
mulation adopted for the rotational broadening (Gray 1 19761 1 
which assumes a rotational profile that is convolved with the flux 
profile. As shown by Conti and Ebbets d 19771 ) this approxima- 
tion may produce a small overstimation of the derived rotational 
velocity of fast rotators. However, as shown by these authors, the 
effect is negligible below 200 km/s and remains small above that 
value. 

The treatment of macroturbulence is a more delicated issue. 
Several authors have pointed out towards indirect evidences of 
some type of macroturbulent broadening acting in the case of 
OB s tars (v iz. Slettebak [19561 Conti & Ebbets [19771 Howarth 
et al. 119971 ); however, its behavior is still not well understood. 
To a first order we will consider an isotropic macroturbulence 
described by a Gaussian distribution of velocities. Some notes 
on other possible descriptions of the macroturbulent broadening 
are given in Sect. [4] 

Finally, in some cases, noise was added to the synthetic 
spectra by using a Poissonian distribution of noise generated 

With RANDOMN in IDL. 

A first test was performed to study how well does this 
method work when different lines from various elements present 
in the stellar spectrum of OB stars are considered. Our main pur- 
pose here is to see whether the Stark broadening affects the use 
of H and He lines as v sin i indicators. This is particularly impor- 
tant when comparing metal and He lines, as the latter have been 
sometimes proposed for v sin i measurements in early type stars 
using the FWHMmethod, particularly when the metal lines were 
faint and shallow, either because of the very early spectral type 
or the large rotational velocity. 

To graphically compare results from different lines, a crA 
baseline was used, as the position of the first zero associated with 
the v sin i depends on A (this is equivalent to a 1 /v scale, see Eq. 
[TJ. Fig. [T] shows the result of this study. The Fourier technique 
has been applied to several Fastwind synthetic lines from dif- 
ferent elements. As it can be seen, the agreement between the 
various lines is excellent (even when the lines He 1/14387 or H^ 
are used) which indicates that using the Fourier method even the 
strongly Stark broadened Balmer lines can be used, in principle, 
as v sin / diagnostics. Note that the power of the Fourier trans- 
form of the He i and H lines concentrates at low frequencies, 
diminishing more rapidly than for the other lines when larger 
frequencies are considered. This is an effect produced by the 
presence of non-rotational broadening mechanisms affecting the 
line profile (i.e. Stark broadening). Similarly to the case of lines 
affected by macroturbulence, it will have consequences in noisy 
spectra. 
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Fig. 1. Four different synthetic line profiles generated with 
Fastwind were convolved with a vsin/ of 60 kms -1 . The dif- 
ferent lines show the ft of the lines on a crA baseline. The first 
zero fits perfectly to its theoretical value for all the lines. The 
ft can detect the vsin / feature even for the He 1/14387 and H^ 
lines. Note that the ft of the H line reaches the value of 1 for 
a- A <200 

How are these results affected by the noise present in real 
spectra? As mentioned by Smith & Gray ( 1976b ), random noise 
produces a "white noise" in the Fourier space (i.e. statistically 
constant with Fourier frequency). Therefore at higher frequen- 
cies, where the power of the signal transform drops, the "white 
noise" dominates and hides the characteristics of the signal 
transform. This effect is illustrated in Fig. [2] A strong synthetic 
O ii line was convolved with two sets of rotational and macro- 
turbulent broadenings: (v sin /, v m ) = (50, 0) and (50, 50) km s . 
Additionally, a Poissonian distribution of noise, corresponding 
to a SNR = 200, was added to both synthetic lines. Fig. [2] com- 
pares the Fourier transforms corresponding to the pure and noisy 
synthetic profiles. For these cases, the white noise level is located 
at ~ 10~ 2 . When the profile is only affected by a rotational broad- 
ening the signal corresponding to the first zero is quite above the 
noise level. This first zero can be perfectly identified, and hence 
the v sin / can be accurately measured. The effect of the noise is 
more dramatic when the line is also broadened by macroturbu- 
lence; in this case, most of the power occurs at low frequencies in 
the Fourier domain and the noise level is above the signal trans- 
form even for the first rotational zero. Consequently, the unam- 
biguous identification of the first zero in the Fourier transform 
becomes difficult or impossible, depending on the combination 
of v sin /, v m and SNR. 

A similar behavior will occur for the H , He i and He n lines. 
Although, in principle, it is possible to use these lines to derive 
the v sin / (as it has been show above), since these lines are broad- 
ened by mechanisms other than rotation, the signatures corre- 
sponding to this rotational broadening may be hidden below the 
noise level due to the lowering of the low frequencies power 
caused by Stark broadening. 

In order to obtain an estimation of the feasibility of the ft 
method when different quality (SNR) spectra are considered we 
constructed a grid of synthetic Si in A 4552 lines convolved with 
various values of v sin / and v m , and degraded to different SNR 
values. The ft method was then used to determine the v sin / for 
the different cases (see results summarized in TableQ]). Although 
we do not expect the noise systematically affecting the position 
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Fig. 2. Effect of the noise on the Fourier transform of a line pro- 
file. A Poissonian distribution of noise generated with idl, and 
corresponding to an SNR = 200, was added to a synthetic On 
line, broadened to a v sin/ = 50 kms~' (up), or to a combina- 
tion of v sin i = 50 km s and v m = 50 km s _1 (bottom). The fig- 
ure compares the transforms of the line with and without noise 
(dashed and solid lines, respectively). 

of the zeroes, this academic example should allow us to better 
understand the applicability limits of the ft method in real spec- 
tra, affected by noise. 



combined with macro turbulence. Although this effect is not 
very important when the macroturbulence value is small, 
deriving correct values even for SNR < 100, the effect may be 
more critical when the amount of macroturbulent broadening 
dominates over the rotational broadening. Note that for the 
case of v sin/ = 60 kms -1 and v m = 80 kms -1 a SNR > 300 
is required to be affected by an error of less than ~ 10% in 
the determination of the v sin i. Finally, Table Q] shows how 
critical the effect of noise can be when the v sin / is low and the 
amount of macroturbulence is large (f.e. v sin/ = 30 kms -1 and 
v m = 80 km s ), a case that we may expect to find in supergiants. 

The final test relates to the effect of microturbulence, which 
may introduce additional zeroes as commented before. This is il- 
lustrated in Fig. [3] For this test a rotational velocity of 60 km s _1 
has been adopted. We see that the Fourier transform of the strong 
O H/14641 line presents zeroes that are different when microtur- 
bulent velocities of 1, 9 and 20 kms _I are introduced. With the 
only information from this line, it would be easy to get a wrong 
value for the projected rotational velocity when using only the 
information from the zeroes (particularly if we only use the first 
one). However, we also see in the figure that we can avoid this er- 
ror by using the information from weaker lines, like O n /14956. 
This line is not significantly affected by the microturbulence, and 
thus only the zeroes due to rotation are simultaneously present 
in both transforms. 
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Fig. 3. Effect of the microturbulence over the determination of 
the v sin i in the Fourier space. 
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Table 1. Formal study of the effect of noise on the v sin i determination 
when a stellar line (Si in /14452) is also affected by macroturbulence. 
We limit ourselves to low values because at higher v sin i the problem is 
much less significant. 



In our formal study, we found that the effect of decreasing 
the SNR is to produce somewhat larger vsin/ values when 



3. Applying the Fourier method to real OB stellar 
spectra 

In the previous section it has been shown that the ft technique 
allows us a straightforward determination of projected rotational 
velocities. When dealing with real spectra one has to remind two 
important properties of the FFT algorithm: 

1. The sampling of the discrete Fourier transform (Act) in- 
versely depends on the product of the number of points into 
which the line profile is sampled and the spectral disper- 
sion: Act- (NAA)~ l . Therefore, the better the sampling of 
the Fourier transform of the line profile, the more accurate 
the position of the zeroes associated to v sin / can be deter- 
mined. 
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Fig. 4. Minimum v sin i that can be determined at different wave- 
lengths for various spectral dispersions as a consequence of the 
combination of formula[T]with the Nyquist frequency 

2. The FFT of a sequence of N real numbers is symmet- 
ric and hence the maximum cr value that is computed is 
cr max = (N/2)(NAAy l = (2AA)' 1 , also known as Nyquist fre- 
quency. From formula[T]and this computational property of 
the FFT it is found that the spectral dispersion (AA) of the 
stellar spectrum imposes a limit in the lowest vsirn that can 
be derived, given by 1.320 c LAI A. Fig.|4]illustrates the min- 
imum v sin i that can be measured for various characteristic 
AA values. 

Therefore, if the spectrum has enough resolution it will be 
possible to determine the projected rotational velocity of the star 
independently of any other broadening mechanism affecting the 
line profiles. The applicability of the ft method in real spectra is 
mainly limited by the spectral dispersion (for low v sin i cases) 
and the SNR (mainly affecting those stellar lines broadened by 
a non-rotational broadening mechanism, i.e. macroturbulence 
and/or Stark broadening). 

Having these ideas in mind, a sample of stellar spectra of OB 
stars, obtained during various observing runs, was restored from 
our archive to determine their projected rotational velocities by 
means of the ft method. The derived values were compared with 
those resulting from the use of the fwhm method. This allowed 
us to test the validity of the ft method by studying those cases in 
which stellar lines are mainly broadened by rotation, as well as 
the range of applicability of the fwhm method. Additionally, the 
possible presence of macroturbulence could be studied. 

3. 1. The fwhm method 

The fwhm method requires the use of an unbroadened line tem- 
plate that is convolved with rotational profiles corresponding 
to various vsin/ values. Several authors have considered dif- 
ferent possibilities for the unbroadened profiles, from sharp- 
lined stars observed with the same configuration than the un- 
knowns (viz. Slettebak |1956l who used 10 Lac; Penny [l996l who 
used HD 34078; and Howarth et al. 119971 who used tSco), to 
synthetic lines from stellar atmosphere code calculations (viz. 
Slettebak et al. 119751 Conti & Ebbets 119771 ), or even Gaussian 
profiles with the same ew than the observed line (viz. Herrero et 



al. 1 1992| l. By determining the fwhm from a Gaussian fit to the ro- 
tationally broadened line templates a fwhm - v sin i calibration is 
then generatecQ. In this study we considered the same approach 
than Herrero et al. ( 11992| l. 

3.2. The observations 

We have selected a sample of Galactic OB stars from observa- 
tions carried out with the Isaac Newton 2.5m and the William 
Herschel 4.2m telescopes at the Observatorio de El Roque de los 
Muchachos during various observational campaigns between 
1989 and 2003. The Intermediate Dispersion Spectrograph 
(ids) and the medium-resolution isis spectrograph attached to 
these telescopes were used. The characteristics of the different 
observing runs are summarized in Table [2] 

The spectra were selected to have SNR > 200 and spectral 
dispersion < 0.4 A/pix (which implies a minimum vsin/ of 
~30-40 kms" 1 , see Fig.|4]l. Only for late O-type and early B- 
type dwarfs some spectra were selected to have a SNR slightly 
below 200, as these targets are not expected to be affected by 
macroturbulence and hence the effect of noise is less critical 
when applying the ft method to metal lines. The sample of stars 
ranges in spectral types from 03 to B2, including dwarfs, giants 
and supergiants (see Tables [3]to[9]and Table ITOb. 

The characteristics of the optical spectra of OB stars do not 
make it possible to use the same set of stellar lines to determine 
the projected rotational velocities for the whole sample of OB 
stars. For example, while in late O-type and early B-type stars 
many metal lines from different elements are present in the stel- 
lar spectra, along with H and He i lines, for the earliest O-type 
stars the spectrum is mainly dominated by H and He n lines, and 
only a few faint metal lines can be found. Therefore we decided 
to perform the study separately following the classification pre- 
sented in Sects !3.3l to 13.61 

3.3. Early B and late O dwarfs: an easy task 

Early B and late O dwarfs have many strong metal lines, along 
with some He i lines, useful to infer the stellar v sin i. We do not 
expect important broadening mechanisms other than rotation 
affecting the metal lines; however, the He lines are additionally 
broadened by the Stark broadening. Tables [3] and [4] show the 
vsin / derived from several Hei, Sim, Cn, Nn, and On lines 
present in the spectra of dwarf stars with spectral types between 
Bl and 09. Results obtained through both the ft and fwhm 
methods are presented. A comparison of v sin i results is shown 
in Figures [5]to[7] (dwarfs are presented as black diamonds). 

For dwarfs, three main conclusions arise from this comparison: 

- He i lines give results consistent with metal lines when ana- 
lyzed with the ft method (Fig. [5] left), but tend to produce 
larger values for the projected rotational velocities when us- 
ing the fwhm method (10-20%, see Figure [5] right). For 
very large rotational velocities, when no metals can be used, 
He i lines give slightly larger vsin / values when using the 
fwhm method than when using the ft, although the effect is 
smaller than for lower v sin i (see Fig. |6j. 

2 Note that there is also the possibility of generating a fwhm - v sin i 
calibration based in previous v sin i determinations (viz. Abt et al. 2002 
Strom et al. |2005t . 
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5500 


5900 A 






ISIS-R1200B 


0.23 


9000 


6500 


6700 A 


CASPEC 




CASPEC® ES03.6 


0.10 


18000 


4050 


5000 A 



Table 2. Observations used for the present work. First column identifies the campaign by telescope and year; second column gives the observation 
date. The different instrument configurations used in each campaign, along with their spectral dispersion, approximate resolution, and observed 
spectral range, are specified in columns third to sixth, respectively. 
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Fig. 6. Comparison of vsini results derived from Hei lines 
through the ft and fwhm methods for v sin i > 180 kms -1 . In this 
range of projected rotational velocities, no metal lines are acces- 
sible for the v sin i determination. 



- From Fig.|7]it can be seen that for 09 - B0 dwarfs both meth- 
ods, the ft and fwhm, give similar results when metal lines 
are used. 

- The shape of the Fourier transform of the metal lines is very 
close to that resulting from a pure rotational profile, confirm- 
ing the absence of important broadening mechanisms other 
than rotation. This fact allows an accurate determination of 
the v sin i of late-O and early-B dwarfs even for SNR ~ 120. 

There are two cases (CygOB2-145, and rSco) in which the 
measurement of v sin i is limited by the spectral resolution. The 
v sin i value resulting from the ft analysis is very close to the 
limit imposed by the Nyquist frequency (i.e. 1.320 cAA/A). In 
addition, the v sin i derived from the fwhm method is somewhat 
below the ft value. A larger spectral dispersion (AA) is needed 
for these stars to correctly derive the vsin/ through the ft 
method. Moreover, even with better resolution spectra, the 
broadening produced by microturbulence may be comparable 
to that one corresponding to the rotational velocity, and hence 
zeroes associated with microturbulence may be found at fre- 
quencies somewhat smaller than cri in the Fourier space. These 
stars may be good examples to study the effects illustrated in 
Fig.S 
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Fig. 7. Comparison of vsin / results derived from metal lines 
through the ft and fwhm methods. Dotted lines represent dif- 
ferences between both determinations of 20 and 40 km s -1 . 



Although we are not going to treat it here in detail, we would 
like to mention that some studies have shown that in the cases 
of large v sin i, the effect of differential rotation may alter the 
position of zeros in the Fourier transform. We refer the reader 
to Reiners & Royer (2004) for a detailed study of this effects in 
the case of A-type stars. However, note that Reiners & Royer 
only find evidence of differential rotation in 4 out of 78 stars of 
spectral types later than ours, and therefore we don't expect this 
effect to be important in our sample. 

3.4. Early B and late O supergiants: the presence of extra 
broadening 

The logical extension of the previous analyses points towards 
early B and late O supergiants since these stars also have many 
narrow metal lines. Results for a sample of supergiants with 
spectral types ranging from 09 1 to B2 1 are presented in Table 
[5] and in Figures to [7] the results from both the ft and fwhm 
methods are again compared (supergiants are shown as red 
triangles). 

In this case of supergiants we find: 

- The ft method gives again consistent results for both metal 
and He i lines (Fig. [5]). 
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Fig. 5. Comparison of v sin i results derived from metal and He i lines through the ft and fwhm methods. Only those cases in which 
v sin i could be derived through metal lines are presented. Dotted lines in the fwhm diagram represents a difference between both 
determinations of 10 and 20 km s _1 . 



Star 


u> 1 Sco 


HD 217657 


HD 37020 


HD 37023 


HD 37042 


rSco 


CygOB2-145 


HD214680 


SpT 


B1V 


B0.5V 


B0.5V 


B0.5V 


B0.5V 


B0.2V 


09.5V 


09V 


Obs. run 


WHT03 


INT92 


INT02 


INT02 


INT02 


CASPEC 


WHT03 


WHT03 


SNR 


350 


120 


220 


290 


350 


250 


120 




Hei!4713 


102(115) 


31 (38) 


53 (58) 


50 (69) 


31 (46) 


14 () 


35 (35) 


36 (48) 


He ii4922 


104 (178) 


b 


b 


b 


b 


b 


b 


b 


Heil5015 


101 (114) 


29 (35) 


60 (63) 


48 (77) 


30 (45) 


15 () 


28 (34) 


35 (50) 


He i 1 5047 


102 (125) 


32 (36) 


52 (55) 


52 (70) 


32 (39) 


18 () 


28 (30) 


33 (46) 


He i ,15875 


104 (119) 












48 (59) 


47 (64 ) 


He i 1 6678 


105 (120) 


29 (49) 


62 (70) 


55 (80) 


35 (53) 




34 (50) 


36 (58 ) 


Si m ,14552 


101 (110) 


36 (21) 


56 (56) 


52 (72) 


32 (36) 


10 (13) 


40 (27) 


34 (40) 


Si ml 4567 


104(110) 


31 (33) 


49 (58) 


50 (69) 


33 (37) 


12(14) 


41 (20) 


38 (46) 


Si mi 4574 


103 (114) 


39 (35) 


53 (54) 


54 (70) 


31 (33) 


11(13) 


39 (24) 


32 (43) 


Cn/14267 


97 (100) 


41 (36) 


49 (53) 


54 (75) 


32 (37 ) 


21 (22) 


36 (38 ) 


40(44) 


N ii A 3995 


103 (102) 


44 (28 ) 












45 (49 ) 


N ii A 4253 


96 (99) 


43 (44) 


56 (60) 


55 (66) 


36 (44) 


b 


43 (36) 


38 (50) 


Nn 14601 


104(110) 


r 


53 (50) 


54 (62) 


b? 


b 


d 


35 (42) 


On/14317 


b 


36 (35) 


56 (54) 


51 (55) 


37 (38) 


14(12) 


44 (31) 


d 


On!4319 


b 


36 (36) 


58 (54) 


47 (53) 


38 (35) 


12 (14) 


39 (25) 


39 (42) 


On/14414 


b 


32 (25) 


49 (52) 


48 (50) 


40 (36) 


14(13) 


35 (29 ) 


d 


On/14416 


b 


30 (29) 


52 (49) 


47 (61) 


33 (36) 


10 (14) 


r 


d 


On A 4590 


108 (109) 


32 (28) 


54 (55) 


48 (63) 


32 (35) 


12(13) 


40 (48) 


b 


On 14595 


97 (105) 


35 (36) 


56 (54) 


48 (52) 


35 (37) 


12 (16) 


42 (24) 


39 (44) 


On 14661 


102 (99) 


28 (24) 


55 (51) 


51 (55) 


30 (33) 


14(13) 


33 (18) 


37 (37) 


On!4699 


b 


29 (35) 


51 (56) 


48 (54) 


35 (38) 


b 


33 (34 ) 


37 (44) 


On!4705 


b 


30(21) 


48 (52) 


44 (52) 


33 (33) 


14 (14) 


35 (28 ) 


34 (37) 


On 14906 


100 (100) 


d 


57 (52) 


46 (51) 


40 (37) 


16 (12) 


d 


d 


On!4943 


b 


29 (26) 


53 (57) 


b 


34 (30) 


12(13) 


29 (21) 


d 



i'sini(FT-He) 103±2 30±2 57 ± 5 51 ±3 33±4 16±2 35 ± 8 37±6 

v sin i (FW-He) 132± 24 40 ±6 62 ± 7 74±5 46±6 ± 42± 12 53±8 

vsini(FT-M) 101 ±4 34±5 53 ± 3 50±3 34±3 13 ± 3 38±4 37±3 

v sin i (FW-M) 105 ± 5 31 ±6 54±3 60±8 36±4 14±3 28±8 43±4 



Table 3. Projected rotational velocities (vsini in kms~') derived from several Hei, Sim, Cn, Nn, and On lines present in the spectra of dwarf 
stars with spectral types between Bl and 09. Values in brackets represent fwhm results. Mean values for He i and metal lines are also presented, 
with errors showing the standard deviation of line values, (b) blended line, (r) noise is important, (d) faint line, (-) line not present in the observed 
spectrum (either because that line does not appear for this SpT, or because it is out of the observed spectral range). 



- For metal lines, the v sin i derived through the fwhm method 
is systematically larger than that resulting from the ft 
method (Fig.[7]i. This result confirms the presence of an ex- 
tra broadening mechanism in early B-type supergiants, as 
shown by Dufton et al. (2006a), and extends the results to- 
wards late O-type supergiants. 

- An additional indication of the presence of an extra broad- 
ening is clearly shown by the shape of the Fourier transform 



of metal lines (see e.g. Fig. [8J. In this case, the power of the 
Fourier transform concentrates at higher frequencies when 
compared to the Fourier transform of a pure rotational line 
profile with the same EW. 

In Sect. 12. H it was shown that the effect of the white noise in 
the Fourier space can be important when an extra-broadening is 
present. If the SNR of the spectrum is not good enough, the zeros 
associated to the v sin ; may be hidden below the white noise, and 
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Star 


HD i 1 Uo 1 




L.ygUriZ-ZJ 


n 11 1 Q /IfiOO 

rSU-1 J4VZ0 


un 1 Anna 
HD 14V / J I 


un 1H(\A 1 


SpT 


B 1 V 


B0.5V 


09.5 V 


09.5 Vn 


(~\C\ 

Uv.jvnn 


09V 


Obs. run 


UN i uz 


rVTT'OO 

1JN 1 ov 


W ±11 Uj 


Vvrll Uj 


Wnl Uj 


UN 1 UZ 


SNR 


250 


1 00 


1 90 


230 


500 


230 


Heil4387 


213 (235) 


109(146) 


158 (185) 


404 (433) 


360 (440) 


132 (180) 


Heil4471 


213 (270) 


91 


159 (212) 


b 


377 (433) 


141 (200) 


Heil4713 


220 (224) 


99(114) 


167 (160) 


d 


b 


127 (145) 


He i 1 4922 


218 (250) 


107 (159) 


177 


405 (416) 


389 (402) 


125 (170) 


Hei/15015 


216 (206) 


- 


153 (162) 


401 (415) 


410(411) 


135 (138) 


He i 1 5047 


(212) 




d 


d 


395 (420) 


130 (166) 


He i A 5875 






154(166) 


429 (421) 


385 (387) 




He i -1 6678 


220 (219) 




164(176) 


b 


404 (407) 


133 (158) 


Sim ,14552 


d 


114(116) 


d 


d 


d 


d 


Si mi 4567 


d 


118(127) 


d 


d 


d 


d 


Si ml 4574 


d 


103 (120) 


d 


d 


d 


d 


Si iv 14089 


d 


b 


153 (157) 


b 


b 


133 (140) 



v sin i (FT-He) 217 ±3 102 ±8 162 ±8 410 ± 13 389 ±17 132 ±5 

vsini(FW-He) 231 ± 23 140±23 177±20 420±8 414±18 165 ± 21 

vsini(FT-M) - 112±8 153 - - 133 

vsini(FW-M) - 121 ±6 157 - - 140 



Table 4. Same as in Table[3]but for stars with v sin i values above 110 km s 1 . For these cases, the high v sin i makes metal lines to be faint and 
blended, and hence not suitable to measure projected rotational velocities. Except for a few metal lines, only He i lines can be used. 



may be mistaken with artificial zeros associated with this white 
noise (see Fig.|2|. A low SNR may result in larger v sin i derived 
values. 

In Fig. [8] the line Sim A 4552 in the isis@wht spectrum of 
HD 190603 (B1.5Ia + ) is used as illustrative example of the appli- 
cation of the ft method to an early B-type Supergiant. The spec- 
trum has a SNR ~ 570 (see Table[3]l. The position of the first zero 
in the ft of the observed line corresponds to a v sin i = 50 km s -1 . 
The shape of the ft of a pure rotational (+ instrumental!) profile 
corresponding to that value does not fit correctly the observed 
case for low frequencies. An extra macroturbulent broadening, 
v m , is needed to produce a satisfactory fit in the Fourier space. 
Note that in this case the power of the ft is below the white noise 
level of the observed line, hence indicating that the ft of the ob- 
served profile may be affected by the white noise, even for those 
frequencies where the first zero associated with the rotation is 
located. The actual v sin i value may be somewhat smaller (and 
consequently, v m larger). 

3.5. Early and mid-0 dwarfs and supergiants 

When moving to earlier spectral types (mid and early O-type) 
most of the strong metal lines present in the optical spectra 
of early B-type stars become very faint or even do not appear 
anymore. A few lines associated with high ionized states of 
C, N, O, and Si can be found in mid O-type stars, and a few 
metal lines in emission are present in the spectra of the earliest 
spectral types. Contrarily, the spectra of mid O-type stars is 
dominated by strong Hei and Hen lines. And finally, even Hei 
lines become weak for the hottest stars. Therefore, in the optical 
spectra of those stars, only a few He lines are available for the 
v sin i determination. 

In previous sections, it has been shown that He i lines result 
in somewhat larger v sin i values than metal lines when measured 
through the fwhm method. Contrarily, the ft method produces an 
excellent agreement. Therefore, application of the ft method to 
the He i lines in O-type stars is very promising to derive reliable 
vsirn values. 

It is important to note that SNR imposes an important limita- 
tion when the ft method is applied to He i lines. Since the Stark 



X \ 


1 1 i 






\ N 








\ v 

\ 


V 






White noise level > 


\\ 

\ ~ 

'" Vi 


AA 






'" u 














HD1 90603 (SHII4552) 


t V ' 






FT observed line 


v. 

1 






FT pure rot profile 


1 






FT rot + macro profile ^ 


t 

. . 1 







0.1 1.0 
o (A -1 ) 

Fig. 8. Fourier transform of the Si III A 4552 line in the isis@wht 
spectrum of HD190603 (B1.5Ia + ). The ft of the observed line 
(solid black) is compared with two synthetic profiles: (dashed 
red) pure rotational profile with v sin; = 50 kms -1 ; (dashed- 
dotted blue) rotational + macroturbulent profile with (vsin/, 
v m ) = (50, 45) kms -1 . Black arrows show the position of the 
zeroes associated with vsin / values of 50 and 68 kms -1 , corre- 
sponding to the values derived through the ft and fwhm methods, 
respectively. 



broadening mechanism has an important influence on the Hei 
lines, high SNR spectra are necessary to be sure that white noise 
in the Fourier space is not affecting the identification of the first 
rotational zero (as shown in Sect. 12.11 ). The effect of noise may 
be even more dramatic if macroturbulence is also affecting the 
He i line profiles. However, note that macroturbulent broaden- 
ing is expected to be more important in giants and supergiants, 
while in these stars Stark effect contributes less to the broaden- 
ing of the Hei lines (i.e. Stark broadening scales roughly as p 2 ). 
Note also that, although the spectra of early O-type stars show 
quite large He n lines, the effect of Stark broadening is larger for 
these lines than for He i lines. Therefore, considering the com- 
ment above the use of these lines will require spectra with very 
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Target 


HD 210809 


CygOB2-10 


HD 209975 


HD 167264 


HD 37128 


HD 38771 


HD2905 


SpT 


091 


09.51 


09.5Iab 


BOIa 


BOIa 


B0.5Ia 


BC0.7Ia 


Obs. run 


INT89 


WHT03 


INT89 


WHTB 


WHTB 


WHTB 


WHTB 


SNR 


200 


250 


200 


400 


330 


400 


300 


C ii 4267 




r 


- 


73 (88) 


75 (94) 


67 (80) 


69 (78) 


Nn 3995 




d 


- 


67 (72) 


64? (89) 


58 (80) 


58 (66) 


Nn 4254 




r 


- 


75 (94) 


69 (92) 


57 (89) 


64 (83) 


Nn 4601 




- 


- 


- 


81 (91) 


59 (89) 


72 (81) 


Nn 4614 




- 


- 


- 


- 


- 


- 


Nn 4621 




- 


- 


- 


- 


- 


- 


N ii 4630 




79 (80) 


69 (71) 


65 (89) 




55 (94) 


62 (78) 


Si ni 4552 




63 (93) 


r 


73 (100) 


83 (98) 


53 (98) 


64 (87) 


Si m 4567 




76(101) 


r 


76(100) 


83 (98) 


52 (96) 


55 (87) 


Si hi 4574 




81 (99) 


r 


77 (103) 


r (77) 


r (94) 


52 (81) 


c: ... A Aon 
Si IV 4089 


1 AC / I ICi 

105 (1z5) 


78 (116) 


59 (110) 


72 (114) 


/9 (1 12) 


73 (99) 


r (86) 


Si IV 41 16 


110(115) 


72 (93) 


70 (108) 


(101) 


65 (98) 


51 (80) 


67 (70) 


On 4590 




76 (106) 


- 


75 (92) 


81 (94) 


51 (86) 


66 (81) 


On 4595 




r 


- 


52 (82) 


81 (94) 


52 (85) 


58 (82) 


On 4661 




61 (76) 


d 


73 (82) 


74 (96) 


53 (86) 


55 (81) 


He 1 4143 






70 (115) 


68 n 001 


78 (941 


68 (96) 


69 (83) 


He 1 4387 


108 (140) 


70 (102) 


80(132) 


75 (107) 


79 (104) 


71 (99) 


65 (99) 


He 1 4471 


113 (150) 


76(113) 




73 (120) 


74 (127) 


b 


b 


He 1 47 1 3 


111 / 1 T7\ 
1 1 1 (12/) 


61 (93) 


75 (117) 


65 (88) 


67 (85) 


55 (78) 


59 (81) 


He 1 4922 


110 (130) 


76 (1 14) 












He 1 50 15 




60 (98) 












He 1 5047 




61 (85) 












He 1 5875 




r 












He i 6678 




r 












v sin i (FT-He) 


111 ±2 


67 ±8 


78±4 


71 ±5 


73 ±6 


63 ±11 


62 ±3 


vsini'(FW-He) 


137 ± 10 


100 ± 11 


125 ±11 


105 ±16 


105 ±21 


89 ±15 


90 ±13 


vsini (FT-M) 


108 ±4 


73 ±7 


66±6 


73 ±4 


76 ±7 


55 ±5 


62±6 


vsini(FW-M) 


120 ±7 


95 ±13 


103 ± 10 


93 ± 11 


94±8 


88±7 


79 ±6 




Target 


CygOB2-2 


HD 13854 


HD 190603 


HD 14956 


HD 14143 


HD 14818 


HD 194279 


SpT 


B1I 


Bllab 


B1.5Ia+ 


B1.5Ia 


B2Ia 


B2Ia 


B21a 


Obs. run 


WHT03 


WHTB 


WHTB 


WHTB 


WHTB 


WHTB 


WHTB 


SNR 




250 


570 


450 


450 


450 


450 


C ii 4267 


43 (50) 


54 (78) 


49 (54) 


(69) 


42 (66) 


68 (76) 


55 (69) 


Nn 3995 




48 (79) 


51 (64) 


46? (72) 


63? (70) 


48 (69) 


54 (74) 


Nn 4254 


59 (64) 


57 (77) 


52 (62) 


56 (67) 


51 (64) 


49 (63) 


55 (73) 


Nn 4601 


68 (68) 


61? (81) 


53 (60) 


55 (77) 


54 (66) 


52 (73) 


52 (74) 


Nn 4614 




60 (68) 


44 (53) 


53? (68) 


56 (59) 


50 (52) 


53 (67) 


Nn 4621 




b(73) 


49 (54) 


52 (69) 


52 (67) 


50 (68) 


52 (66) 


Nn 4630 




52 (84) 


52 (63) 


54? (76) 


53 (73) 


54 (70) 


59 (74) 


Sim 4552 


44 (58) 


r(85) 


50 (68) 


55 (76) 


48 (70) 


47 (74) 


49 (75) 


Sim 4567 


48 (61) 


45 (83) 


44 (68) 


55 (75) 


54 (72) 


54 (74) 


52 (71) 


Sim 4574 


45 (55) 


54 (77) 


52 (63) 


61 (68) 


(66) 


60 (68) 


52 (68) 


Si iv 4089 




r(73) 


b 


d 


d 


d 


d 


Siiv4116 


58 (56) 


65 (73) 


49 (59) 


d 


d 


d 


d 


On 4590 


42 (67) 


r(74) 


57 (59) 


61 (69) 


44 (66) 


57 (64) 


57 (67) 


On 4595 


44 (63) 


r(67) 


55 (56) 


43 (69) 


46 (54) 


53 (67) 


56 (71) 


On 4661 


47 (41) 


r(76) 


52 (58) 


48 (69) 


47 (59) 


54 (68) 


48 (67) 


Hei4143 


b 


62(86) 


54 (67) 


48 (80) 


53 (71) 


57 (81) 


56 (76) 


He i 4387 


50 (92) 


48(90) 


46 (74) 


48 (87) 


53 (81) 


61 (85) 


59 (82) 


He i 4471 


b 


62(110) 


59 (88) 


73 (98) 


67 (98) 


74(103) 


66 (95) 


Hei4713 


47 (63) 


56 (77) 


50 (62) 


47 (69) 


52 (66) 


64 (72) 


48 (65) 


He i 4922 


b 














Hei5015 


49 (67) 














He 1 5047 


49 (58) 














He 1 5875 


59 (82) 














He i 6678 


50 (79) 














V sin i (FT-He) 


51 ± 5 


55 ±7 


52±7 


48±1 


53 ±1 


63±2 


54±8 


vsini(FW-He) 


74 ± 13 


92 ±17 


75 ±13 


85 ±15 


82 ± 16 


87 ± 16 


87 ±15 


v sin i (FT-M) 


49±9 


55 ±6 


51 ±4 


53 ±5 


51 ± 6 


54±6 


53 ±3 


ii sin i (FW-M) 


58 ±8 


77 ±6 


60±5 


71 ±4 


66±6 


69±6 


71 ±3 



Table 5. Same as Table[3]but for Early-B and late-0 supergiants 



high SNR. We did not achieve enough SNR in any of the spectra 
of our sample of early O-type stars, and hence He n could not be 
used. 

Tables [6] and [7] summarize the line-to-line vsin/ values 
resulting from the fwhm and ft analyses of the sample of early 
and mid O-type dwarfs and supergiants. Hei lines were mainly 



used; however, for the supergiants, we could also use some 
metal lines (some of them present in the spectral range between 
5000 and 6000 A, only observed in the wht03 campaign). A 
comparison of v sin i results is again shown in Figures [5] to [7] 
(dwarfs and supergiants corresponding to black diamonds and 
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Star HD 12993 CygOB2-16 HD46056 HD 13268 HD46966 HD236894 HD 168137 

SpT 06.5V 07.5V 08Vn ON8V 08V 08V 08.5V 

Obs.run PerOBl WHT03 CAHA INT89 INT89 PerOBl WHT03 
SNR 



He i A 4387 77 (123) 

He i i 4471 80(90) 

Heii4713 87(122) 
He i i 4922 
He i A 5875 

Si iv A 4089 84(90) 

Sirv-14116 r 

O in A 5592 

C iv A 5801 
Civ /1 5811 



209 (220) 345 (352) 
214 (230) 350 (404) 
203 (216) 

191 (217) 359 (371) 
201 (224) 

b 
b 

213 (238) 
b 

224 (235) 



315 (307) 70(100) 
310(320) 54(98) 
57 (72) 

298(301) 54(80) 



148(170) 80(106) 
142(176) 59(110) 
148 (156) 48 (73) 
58 (100) 
53 () 

144 (152) 53 (79) 
b 61 (69) 

58 (80) 

50 (77) 
50 (77) 



V sin i (FT-He) 81 ±5 

v sin / (FW-He) 112 ±9 
v sin i (FT-M) 84 
v sin i (FW-M) 90 



204±9 351 ±7 

221 ±6 376 ±26 
219±8 
237 ±2 



308 ± 9 59 ± 8 

309 ± 10 88 ±15 



146 ±3 61 ±13 

167 ±10 97 ±17 

144 54 ±5 

152 76 ±4 



Table 6. Same as Table[3]but for mid-0 dwarfs 



Star 


CygOB2-7 


CygOB2-8C 


CygOB2-9 


CygOB2-ll 


CygOB2-8A 


HD 210839 


HD 192639 


SpT 


03If 


05If 


05If 


05If+ 


05.51(f) 


06I(n)fp 


07Ib(f) 


Obs. run 


WHT03 


WHT03 


WHT03 


WHT03 


WHT03 


INT92 


INT89 


SNR 


250 


250 


150 


200 


350 


150 


220 


He 1/14387 














95 (121) 


He 1/14471 




180 (211) 


97 (136) 


107 (154) 


87 (136) 


223 (256) 


104(163) 


Heii4713 










84 (102) 




114(125) 


He i A 4922 














115 (130) 


He i A 5875 


94(117) 


170 (205) 


97 (119) 


120(134) 


82 (125) 






Cm A 5695 




(249)" 


104 (120)" 


120 (140)" 


98(118)" 






C iv A 5801 


96 (100)° 


b 


b 


b 


b 






Civ ,15811 


105 (103)° 


180 (200) 


b 


117 (143) 


84 (120) 






O in A 5592 




184 (196) 


113 (135) 


100 (136) 


80(135) 






Si iv A 4089 


89 (112)" 












100(114) 


SiiV/14116 


89 (112)" 


(233)" 


(179)" 


122(147) 


(108)" 




117 (112) 


NnU4634 


(120)° 


b 


b 


b 


b 


b 


b 


N mA 4641 


94 (120)" 


b 


b 


b 


b 


b 


b 


N iv A 4058 


102 (143) 















223 107 ±9 

256 135 ±19 

109 ± 12 
113 ± 1 



Table 7. Same as Table[3]but for early and mid-0 supergiants. " Line in emission. 



v sin ; (FT-He) 94 

v sin i (FW-He) 117 

v sin i (FT-M) 96 ± 7 

v sin i (FW-M) 116 ±14 



175 ±7 97 

208 ±4 128 ±12 

182±3 109±6 

220 ±26 145 ±31 



114 + 9 84±3 
144 ±14 121 ±17 

115 ± 10 87 ±9 
142 ±5 120 ±11 



red triangles, respectively). 

From the comparison of v sin i values obtained for metal lines 
with the ft and fwhm methods in some of the supergiants in our 
sample, we can conclude that macroturbulence is also present 
in these stars, extending the result found in B-type supergiants 
to the O-type supergiants. Unfortunately, the presence of macro- 
turbulence in the early O-type dwarfs could not be investigated 
because there was no metal lines available for the v sin i determi- 
nation. 

3.6. OB giants 

Results for the determination of projected rotational velocities 
in the sample of giant stars are presented in Tables [8] and [9] 
Similar to previous sections, v sin i measurements from both the 
ft and fwhm for several lines, as well as mean values for He i 



and metal lines, are summarized in those Tables. For these stars 
He i lines were mainly used, although in some cases a few Si iv, 
C iv, N ii and O m lines were also available. 

The analysis of the sample of Giant stars points towards 
two curious results: 

- Both low and high v sin i values are found. 

- Macroturbulent effect is either very large or negligible. 

These will be discussed in more detail in the next section. 

4. Discussion 

Fig.|5]shows a comparison of v sin i results obtained for the sam- 
ple of stars when applying the ft and fwhm methods to metal 
and He i lines. This comparison was only possible for stars with 
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Star 


CygOB2-8B 


CygUBz-4 


HD 1 64492A 


T TT~\ 1 C O C n A 


SpT 


06.5111 


u/ni(i) 


07.5III 


08III 


Obs. run 


11 7"T TTA *") 

WH103 


WH 1 03 


IN 102 


WH103 


SNR 


230 


220 


250 


300 


He i A 4387 


d 


d 


45 (89) 


88 (112) 


Heii4471 


74 (130) 


83 (120) 


49 (94) 


79 (129) 


l|„ T } AIM. 

rle I/f+z 13 




Q7 (\ 1 S\ 
oz ^1 1 




77 


He iA 4922 




(128) 


52 (80) 




He i i 5015 


94(115) 


85 (116) 


48 (60) 




He a 5875 


85 (115) 


80(130) 






Civ A 5811 


83 (110) 


70 (136) 






N n/14254 


65 (118) 


74(125) 






in A 5592 


82(118) 


74(125) 






Si iv A 4089 


84(119) 


84(125) 


50 (61) 




SirvA4116 




81 (103) 


45 (56) 


85 (89) 


Siivi4212 






50 (51) 


81 (81) 


Si iv A 4654 






45 (51) 


76 (81) 




v sin (' (FT-He) 


84 ± 10 


83 ±2 


49±3 


81 ±6 


v sin i (FW-He) 


120 ±9 


122 ±6 


77 ±16 


110 ±20 


vsini(FT-M) 


79±9 


77 ±6 


48 ±3 


81 ±5 


v sin I (FW-M) 


116±4 


123 ± 12 


55 ±5 


84 ±5 



Table 8. Same as Table [3] but for OB giants with 1'sini values below 
100 km s-'. 



500 



400 



300 



.E 200 
> 

100 



~i 1 1 1 1 1 1 1 r 



O Dwarfs 

A Supergiants 

□ Giants 



□ 

O 



O 



□ 



o 
o 



A □ 



o 



A 



O 



o 



A 



A 

a loin 



o 



oo 



□ 



03 04 05 06 07 08 09 B0 B1 
SpT 

Fig. 9. Distribution of projected rotational velocities (as deter- 
mined by the ft method) among the various spectral types of the 
stars in the sample. Both v sin /ft (metal) and vsin/p^Hei) are 
used in this figure. Horizontal dashed line represent the resolu- 
tion limit for R = 9000. 



Star 
SpT 
Obs. run 
SNR 


HD24912 
07III(n)((f)) 
INT89 
250 


HD 203064 
07.5III:n((f)) 
INT89 
200 


HD 191423 
09III:n 
INT89 


He i ,14387 
Heii4471 
Heii4713 
He i ,14922 


224 (247) 
211 (247) 

225 (242) 
215 (230) 


317 (338) 
330 (322) 
328 (335) 


468 (470) 
429 (437) 




v sin r (FT-He) 
v sin i (FW-He) 


219±7 
242 ±8 


325 ±7 
332 ±9 


450 ±30 
454 ±23 



Table 9. Same as Table[3]but for OB giants with large v sin ('. 



vsin/ values below 200 kms . For those cases with a pro- 
jected rotational velocity above this value metal lines appeared 
blended or very faint, and hence they were not available to derive 
v sin i. For the studied stars, the ft method produces coherent re- 
sults between metal and He i lines, indicating the possibility of 
accurately deriving projected rotational velocities through this 
method by using Hei lines. This is not the case for the fwhm 
method, whose He i lines v sin / associated values are somewhat 
larger than those derived from metal lines, even for large values 
of v sin i (see also Fig. [5). This result is a consequence of the 
fact that He i lines are affected by Stark broadening, while we 
are considering a Gaussian definition of the intrinsic profiles. It 
indicates that the assumption made by Herrero et al. (1992), that 
this hypothesis can be used when applying the fwhm method to 
He i lines for determining v sin / if this is large, is not completely 
correct. The effect of pressure broadening over He i may be im- 
portant even for large projected rotational velocities, although 
for v sin i > 250 km s _1 the relative error is small. Note that, al- 
though this discrepancy should not be present if the intrinsic pro- 
files used for the He i lines consider the Stark broadening, the de- 
rived fwhm values may require a large computational effort and 
be model dependent , while ft values will not. 

Since the Stark effect over the Hei lines scales roughly 
as p 2 , these lines are expected to be less affected by Stark 



broadening in supergiants. This explains why for some of the 
studied supergiants, a better Hei vs. metal lines agreement is 
found when applying the fwhm method (Fig. [5] right). However, 
caution!, the agreement between Hei and metal lines does not 
necessarily mean that the fwhm method is giving correct values 
for the v sin i. 

Differences between the v sin i values obtained through the 
application of the fwhm and ft methods to metal lines can be 
interpreted as an indicator of the presence of some extra broad- 
ening, traditionally associated to large scale motions. As indi- 
cated above, we will refer to it as macroturbulence. This is il- 
lustrated in Fig. [7] The whole sample of stars whose metal lines 
are available for the v sin i to be measured are considered in this 
figure. Note that a similar comparison with He i lines (i.e. Fig. 
[6]l does not have sense for investigating the presence of macro- 
turbulence using the differences between the ft and fwhm meth- 
ods since these lines are also affected by other non-rotational 
broadenings as shown above (metal lines are also affected by 
other line broadening mechanisms, like thermal or collisional 
ones, but we expect them to be negligible, as shown by the 
agreement found in dwarf analyses). Stars have been separated 
according to luminosity class. While macroturbulence can be 
considered negligible in dwarfs, supergiants are clearly affected 
by this extra-broadening mechanism; 65% of the sample of su- 
pergiants shows vsin/(Fr)< 0.8 v sin /(fwhm). The case of gi- 
ants is more complicated; these objects display a large range of 
v sin z'(ft)/v sin /(fwhm) ratios. 

Supergiants rotate as a group more slowly than dwarfs 
do. Except for two cases among early O-type supergiants 
(CygOB2-8C and HD 218039) all of them have projected 
rotational velocities below or very close to 100 kms -1 . This has 
to be compared with the large range of v sin i found for dwarfs 
and giants, that cover from about 15 up to 450 km s . In Fig. [9] 
we see that the region of low rotational velocities is shared by all 
spectral types and luminosity classes, but that of large rotational 
velocities is dominated by giants and dwarfs of mid and late 
spectral types. This result was already indicated by Slettebak 
d 19561 ) and later on by other authors (viz. Conti & Ebbets 
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119771 Howarth et al. 11997) . These authors also propose indirect 
indications of the possible presence of macroturbulence in OB 
stars, but could not directly confirm the amount of broadening 
that is actually due to rotation and therefore their v sin /values 
are affected by the contribution of macroturbulence. The ft 
method shows that indeed very large rotational velocities are 
the cause for the large line broadening in dwarfs, and allows 
for a quantitative estimation of the actual v sin i in supergiants 
(where effects of macroturbulence are important). First results 
for a more limited sample were presented by Ebbets ( 1979 ). Our 
results also support the bimodality found by Conti and Ebbets 
d 19771 1 and other authors since then in the v sin i distribution of 
main sequence stars. 

Although it is clear that supergiants are affected by macro- 
turbulence, our knowledge of its behavior and the physical 
mechanisms which produce this extra broadening is still very 
poor. In Sect. 12.11 an isotropic Gaussian macroturbulence was 
assumed for the formal tests; however, as shown by Dufton et 
al. (12006al l. one cannot rule out other macroturbulent velocity 
distributions from the comparison with high resolution observed 
line profiles. Physically, a Gaussian-angle-dependent velocity 
distribution similar to that proposed by Gray (1976 ) to model 
the flow of convective cells in the Sun, may also have sense 
in the case of B supergiants. In that case, the large-scale 
radial-tangential movements may be produced when the stellar 
material flowing outwards as a stellar wind is coupled with 
the stellar rotation. Observational tests to investigate the exact 
degree of the anisotropy of the macroturbulence should be 
investigated in the future. In the meantime, rough estimations 
of the macroturbulence are also valuable, once the projected 
rotational velocity has been inferred. Note that the characteristic 
value of the macroturbulent velocity is dependent on the ap- 
proach considered for the macroturbulent velocity distribution 
(see p.e. Dufton et al. l2006al ). 

The quality of our observational data-set (mainly selected to 
show the applicability of the ft method to OB stars) is not good 
enough to proceed to a fine estimation of the macroturbulenc^- 
However, taking into account the discussion above, we have de- 
cided to consider a rough estimation of the amount of extra- 
broadening, given by ©o, where this quantity is calculated as 
follows. A v sin / -fwhm calibration is created by measuring the 
fwhm corresponding to the convolution of an instrumental pro- 
file and various values of v sin i. This calibration is applied to a 
second set of line profiles in which an isotropic Gaussian macro- 
turbulence (with characteristic velocity ©o) is also included. The 
corresponding v sin i fwhm is hence inferred for various v sin i - 
©o combinations, ending with a set of lines in the v sin / fwhm - 
v sin /ft plane (see Fig. ITQb. Once this two quantities are mea- 
sured for a given star, the ©o can be estimated by interpolation 
between these curves. We expect the quantity ©o calculated in 
this way to be roughly proportional to the macroturbulent veloc- 
ity and behave in a similar way when varying the stellar param- 
eters. 

The influence of ©o relative to v sin i seems to increase when 
moving from dwarfs to supergiants, as can be seen in Fig.[TTJ We 
have decided to plot ©o/v sin / vs. SpT, instead of ©o, because for 
large values of v sin i is more difficult to distinguish between the 
different values of ©o, as it is illustrated by the proximity of the 
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Fig. 10. This plot illustrates the method used to derive ©o, a 
rough estimation of the macroturbulence. The various dashed 
curves represent cases with different values of ©o, from 5 to 95 
km s . An instrumental profile (R = 9000) has also been consid- 
ered (see text for explanation). Symbols have the same meaning 
than in previous figures. 




3 We refer to the works by Dufton et al. d20C )6a) and Lefever et 
al. ( 1200 6) as examples of macroturbulent determination through better 
quality observational data. 



03 04 05 06 07 08 09 B0 B1 
SpT 

Fig. 11. Distribution of the ratio ©o/v sin / in the sample of stars 
as a function of spectral types. Only those stars with metal lines 
available to measure vsin/pw and v sin /ft were used (see text 
and Figure [TO]). 



dashed curves when v sin i increases. We consider this ratio more 
meaningful. 

Except for the two giants around 07, the largest values at 
each spectral type are always obtained for supergiants. Only 
two dwarfs have ©o/v sin/ comparable to the one found in 
Supergiants; these are HD 37023, which is known to be binary, 
and HD 168137 (08.5V). 

Combining both results, we conclude that the role of rotation 
should decrease and the role of macroturbulence increase when 
the stars move from dwarfs to supergiants, i.e, when they evolve. 
This can also be seen in Figures Q~2] and [13] where we have rep- 
resented the stars on a HR diagram, with symbol sizes propor- 
tional to their rotational and macroturbulent velocities. To place 
the stars in the HRD we have taken results from various authors 
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Fig. 12. HR diagram where the stars in the studied sample have 
been located with symbol sizes proportional to their rotational 
(v sin /ft) velocities. As in previous figures, black, blue, and red 
symbols represent dwarfs, giants and supergiants, respectively. 

that used state-of-the-art nlte, line-blanketed, unified model at- 
mospheres (Herrero et al. 120021 Repolust et al. 2004 Urbaneja 
120041 Simon-Diaz 120051 Simon-Diaz et al. 120061 Crowther et al. 
2006 ). For those stars in the sample which have not been previ- 
ously analyzed, we use the SpT-T e ff and SpT-L calibrations de- 
fined by Martins et al. (2005 ). Note that for some of the B0.5V 
stars, the stellar parameters were obtained by extrapolating the 
calibration by Martins et al. 

We see that the symbols in Fig.[T2]may have any size close 
to the ZAMS, but they become smaller when we depart from 
it. Beyond an imaginary line parallel to the ZAMS all symbols 
are very small, indicating slow rotation. Although the inclina- 
tion angle is unknown, it is obvious that the effect cannot be 
attributed to a random effect. Clearly, all supergiants have small 
rotational velocities which, in view of the distribution for dwarfs 
and giants, indicates that the stellar rotational velocities decrease 
when the stars evolve. As suggested by other authors (see review 
by Herrero & Lennon 2004 ) this may be related to the loss of an- 
gular momentum caused by stellar winds. 

Recently Dufton et al. (2006a) and Lefever et al. (2006) have 
applied the ft method to two samples of Galactic B-type super- 
giants observed with high spectral resolution. They found larger 
vsin/ values for the earliest spectral types. We put together the 
vsirn results by these authors and our results in Fig. [14] We can 
see that they agree very well, and that our results extend those of 
Dufton et al. and Lefever et al. towards higher temperatures. 

Additionally, these authors could use the good quality of 
their data to estimate the amount of macroturbulent broadening 
present in the stars by means of the fit of the theoretical line pro- 
files to the observed ones. Since we limit ourselves to a crude 
first order approximation as shown above, we have to be very 
careful in interpreting the results concerning our estimation of 
what we have called macroturbulencfl- However, it is interesting 
to see how this estimation varies across the HR diagram. In Fig. 
[T3lwe see that the symbol sizes representing the macroturbulent 
velocity have a much more uniform distribution than those of 

4 Note that even a comparison of the values inferred by these authors 
for the macroturbulence should be handle with care, since Dufton et al. 
are using an isotropic Gaussian definition while Lefever et al. assume a 
radial-tangential approach. 
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Fig. 13. Estimates of projected rotational velocities for our sam- 
ple of Supergiant stars. Results from the analysis of the Galactic 
B-type supergiants by Dufton et al. and Lefever et al. are also 
included. 
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Fig. 14. Same as in Fig. Q~2] but for the quantity ©o (i.e. rough 
estimation of the macroturbulence) 

the rotational velocities, with the lowest values corresponding to 
dwarfs close to the ZAMS. Part of this difference may be related 
to projection effects (that would affect large scale turbulent mo- 
tions much less than rotation). Nevertheless, according to these 
estimations, macroturbulence does not necessarily increase sig- 
nificantly when the stars evolve redwards, but may remain ap- 
proximately constant. However, its role in line broadening rel- 
ative to rotation increases (and perhaps also in other physical 
mechanisms). Note that this broad scenario suggested here (ro- 
tational velocity decreasing when the stars evolve from dwarfs 
to supergiants, macroturbulence being roughly constant but tak- 
ing over a dominant role) has been deduced from the study of a 
very inhomogeneous sample of OB stars. To better understand 
the role that rotation and macroturbulence play in the evolution 
of massive stars, an homogeneous study of stars in clusters of 
different ages and metallicities should be performed (viz. Dufton 
et al. l2006bl Mokiem et al. 120061 

5. Conclusions 

In the past the Fourier method for the determination of projected 
rotational velocities has been extensively used in the analysis 
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of late-type stars, but only marginally applied to the early-type 
case. In this paper we show the strength of this method for its 
application to OB stars. 

The noise and the spectral dispersion of the observed spectra 
are the main limitations to the applicability of the Fourier 
method. Spectral dispersions below 0.4 A/pix are required to 
measure vsin/ values below 30-40 kms _I when using the 
optical spectrum of OB stars. For dwarf stars, not very affected 
by non-rotational broadenings (i.e. macroturbulence), the ft 
method can be successfully applied even for SNR values as 
low as 100. The effect of noise is more important in the case 
of supergiants, where absorption stellar lines are broadened by 
mechanisms other than rotation. If the SNR of the observed 
spectra is not high enough and the amount of macroturbulent 
broadening is of the order of rotational broadening, the applica- 
tion of the ft may result in somewhat larger v sin i values than 
the actual ones; SNR values above 300-400 are required for 
these cases. 

The comparison of v sin i values obtained through the ft and 
fwhm methods in a sample of O and early-B stars shows that 
the fwhm technique must be used with care in the analysis of 
OB giants and supergiants, and when it is applied to Hei lines. 
Contrarily, the ft method appears to be a powerful tool to derive 
reliable projected rotational velocities, and separate the effect 
of rotation from other broadening mechanisms present in these 
stars. 

The analysis of the sample of OB stars, that includes dwarfs, 
giants and supergiants with spectral types earlier than B2, shows 
that while dwarfs display a broad range of projected rotational 
velocities, from less than 30 up to 450 km s _1 , supergiants 
have in general values close to or below 100 km s _1 . Giants 
also display a large range of projected rotational velocities. 
The analysis has also definitely shown that while the effect 
of macroturbulence in OB dwarfs is usually negligible when 
compared with the rotational broadening, the effect of this 
extra-broadening is clearly present in supergiants. 

When examining the behavior of the projected rotational ve- 
locities with the stellar parameters and across the HR diagram, 
we conclude, in agreement with previous researchers, that the ro- 
tational velocity should decrease when the stars evolve. On the 
contrary, extra broadening across the HR diagram, interpreted 
as macroturbulence, may be constant, resulting therefore in an 
increasing importance as compared to rotation when the stars 
evolve. 
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Star 


SpT 


v sin i'(ft) 


v sin i'(fwhm) 


T c ,r(K) 


log L/Lq 


Ref 


CygOB2-7 


03If* 


96 ±7 


1 16 ± 14 


45500 


5.91 


He02 


CygOB2-8C 


051(f) 


182 ±3 


220 ±26 


41000 


5.66 


He02 


CygOB2-9 


05K+ 


109 ±6 


145 ±31 


38610 


5.87 


Marf)5 


CygOB2-ll 


05I(f)+ 


115 ± 10 


142 ±5 


37000 


5.92 


He02 


CygOB2-8A 


05.511 


87 ±9 


124 ±9 


38500 


6.19 


He02 


HD 210839 


06I(n)fp 


223M 


256<«> 


36000 


5.83 


Rep04 


HD 192639 


07Ib(f) 


109 ± 12 


1 13 ± 1 


35000 


5.68 


Rep04 


HD210809 


09Iab 


108 ±4 


120 ±7 


31500 


5.60 


Rep04 


CygOB2-10 


09.5Iab 


71 ±9 


95 ± 12 


29000 


5.77 


He02 


HD 209975 


09.5Iab 


66±6 


103 ± 10 


32000 


5.69 


Rep04 


HD 167264 


BOIa 


73±4 


93 ± 11 








HD 37128 


BOIa 


76±7 


94±8 


27500 


5.69 


Urb 


HD 38771 


B0.5Ia 


55 ±5 


88±7 


26500 


4.94 


Urb 


HD2905 


BC0.7Ia 


62 ±6 


79±6 


21500 


5.52 


Cro06 




D 1 I 

oil 


AQ _i_ Q 


JO ± 


ZoUUU 


4.0J 


rieuz 


HD 13854 


Bllab 


55 ±6 


77 ±6 


21700 


5.92 


Urb 


HD 190603 


B1.5Ia+ 


51±4 


60±5 


18500 


5.57 


Cro06 


HD 14956 


B1.5Ia 


53 ±5 


71±4 


20500 


5.66 


Urb 


HD 14143 


B2Ia 


51 ±6 


66 ±6 


19500 


5.48 


Urb 


HD 14818 


B21a 


54±6 


69±6 


20100 


5.47 


Urb 


HD 194279 


B2Ia 


53 ±3 


71 ± 3 


19000 


5.37 


Cro06 


PvctOR2-8R 


06 MITff) 


79 ± 9 


1 16 ± 4 


37130 


5.54 




HD24912 


07III(n)((f)) 


219±7 (D) 


242 ±8 <fl) 


35000 


5.42 


Rep04 


CygOB2-4 


07111(f) 


77 ±6 


123 ±12 


35500 


5.41 


He02 


HD 164492A 


07.5111 


48 ±3 


55±5 


35020 


5.42 


Mar05 


HD 203064 


07.5III:n((f)) 


325±7< fl) 


332±9 W 


34500 


5.50 


Rep04 


HD 168504 


08III 


81+5 


84±5 


33960 


5.35 


Mar05 


HD 191423 


09III:n 


450±30 ( °) 


454±23< a) 


32500 


5.23 


Rep04 


HD 12993 


06.5V 


84 


90 


37870 


5.23 


Mar05 


CygOB2-16 


07.5V 


219±8 


237 ±2 


35870 


5.05 


Mar05 


HD 46056 


08V((f)) 


351±7 <fl) 


380 ± 30 <a) 


34980 


4.96 


Mar05 


HD 13268 


ON8V 


308±9 (o) 


309 ± 10 <a) 


33000 


5.05 




HD 46966 


08V 


59 ± 8 <») 


88 ± 14<»> 


34880 


4.96 


Mar05 


HD 236894 


08V 


144 


152 


34880 


4.96 


Mar05 


HD 168137 


08.5V 


54±5 


76±4 






Mar05 


HD 37041 


09V 


132±5<"> 


165 ±21 


35000 


4.93 


SD06 


HD 214680 


09V 


37 ±3 


43±4 


36000 


5.02 


SD06 


HD 149757 


09V 


390 ± 17 ( "> 


414 ± 18 (o) 


32880 


4.77 


Mar05 


BD- 134928 


09.5Vn 


410 ± 13<"> 


420±8 < ° ) 


31880 


4.68 


Mar05 


CygOB2-23 


09.5V 


162±8 (o) 


177±20 <fl) 


31880 


4.68 


Mar05 


CygOB2-145 


09.5V 


38±4 


28 ±8 


31880 


4.68 


Mar05 


tSco 


B0.2V 


<13 


14±3 


32000 


4.47 


SD06 


HD 37020 


B0.5V 


53 ± 3 


54 ±3 


30000 


4.45 


SD06 


HD 37023 


B0.5V 


50±3 


60±8 


32000 


4.47 


SD06 


HD 37042 


B0.5V 


34±3 


36±4 


29000 


4.11 


SD06 


HD217657 


B0.5V 


34±5 


31 ±6 








HD 228 199 


B0.5V 


112±8 


121 ±6 








HD 37061 


B1V 


217±3< fl) 


231±23 (a) 


32000 


4.54 


SD05 


to 1 Sco 


B1V 


101±4 


105 ±5 









Table 10. All vsini values correspond to metal lines measurements except those marked with <a) . References for the effective temperatures 
and luminosities: He02 (Herrero et al. 2002), Rep04 (Repolust et al. 2004), Urb04 (Urbaneja, 2004, thesis) , Cro06 (Crowther et al. 2006), SD05 
(Simon-Diaz, 2005, thesis) SD06 (Simon-Diaz et al. 2006). For the stars in the sample which have not been previously analyzed with line-blanketed 
models, we use the SpT - T eff and SpT - L calibrations defined by Martins et al. (2005, Mar05). 



